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Abstract. Recursion formulae are derived for the calculation of two 
centre matrix elements of a radial function in relativistic quantum 
mechanics. The recursions are obtained between not necessarily 
diagonal radial eigensates using arbitrary radial potentials and any radial 
functions. The only restriction is that the potentials have to share 
a common minimum. Among other things, the relations so obtained 
can help in evaluating relativistic corrections to transition probabilities 
between atomic Rydberg states. 
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1. Introduction 

The bridge between quantum calculations and experimental results is made at the level 
of matrix elements of the appropriate operators (Wong and Yeh 1983a,b De Langc and 
Raab 1991, Dobrovolska and Tutik 1999, Quiney et al 1997, Elboudali and Joulakian 
2001, Ilarraza-Lomeli et al 2002). In atomic and molecular physics this usually means 
matrix elements of powers of a radial coordinate between states of the system at hand 
(Blanchard 1974, Brack 1983, Bessis et al 1985, Nana Engo et al 1997, Owono Owono 
et al 2002, West 2001). Matrix elements of more general radial functions are also 
very useful (Avram and Draganescu 1997, Charro and Martin 2002, Nagy and Vegh 
1992). In nonrelativistic quantum mechanics, the importance of hypervirial results, 
and other related techniques, follows from this fact, since the task of calculating matrix 
elements is simplified using such techniques (Morales 1987, Fernandez and Castro 1987, 
1996, De Lange and Raab 1991, Nunez- Yepez et al 1995). These techniques are also 
important for atomic physics in the relativistic realm (Martinez-y-Romero et al 2001). 
This is especially important at present given the precision attained in atomic physics 
experiments using synchroton radiation and other state-of-the-art sources of photons 
(Schippers et al 2002, 2003, Miiller et al 2002, Aguilar et al 2003). 

In this work we exhibit that hypervirial-like formulas and other operator 
techniques suffice to obtain recurrence relations between relativistic matrix elements 
of an arbitrary radial function between states corresponding to two different arbitrary 
potentials. Such relation generalizes other recurrence relations between relativistic 
hydrogenic matrix elements of powers of r and of /3r [/3 is a Dirac matrix, see 
equation (2)] we have recently calculated (Ilarraza-Lomeli et al 2002, Martinez- 
y-Romero et al 2000, 2001, 2002). In subsections 2.1, 2.2 and 2.3 of the paper 
we obtain the four recursion relations between the matrix elements of an arbitrary 
radial function, /(r), with the corresponding matrix elements of its first and second 
derivatives, taken between eigenstates of different potential functions, Vi{r) and V2{r), 
with the corresponding matrix elements of its first and second derivatives. These 
generalized recurrence relations follow from a hypervirial-like result together with 
operator algebra computations. The recursions so obtained may be useful for studying 
radiative transitions in Rydberg atoms, in analysing atomic photorecombination and 
photoionization processes, for calculating relativistic corrections to ionic oscillator 
strengths, or in analysing impact ionization or vibrational transitions in molecules — 
albeit in the last two cases in a crude manner — (Schippers et al 2002, Owono Owono 
et al 2002, Semenov and Cherepkov 2003, Week et al 2001, Bauche-Arnoult et al 
1989). Section 3 contains the conclusions and other final remarks. 

2. Relations for two centre integrals in relativistic quantum mechanics 

General hypervirial results, the virial theorem, and other algebraic and operator 
techniques have been always very useful for calculating matrix elements in quantum 
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mechanics. Nonetheless they have been httle used in relativistic situations, see, 
however, (Brack 1983, Lucha and Schoberl 1990). We want to show in this paper that 
they have similar importance for obtaining recursion relations between matrix elements 
of radial functions between eigenstates of different radial potentials in relativistic 
quantum mechanics. Thus, let us consider two radial Dirac Hamiltonians with two 
possibly different radial potentials (each behaving as the temporal component of a 
4- vector) Vi(r) and V2(r). We further assume that these potentials have the same 
equilibrium position which, furthermore, is assumed coincident with the origin of 
coordinates. That is, the recurrence relations correspond to the so-called unshifted 
case. The potentials could be regarded as describing the electronic configurations in 
processes involving electronic rearrangements in atomic transitions. For example, any 
transition to an autoionizing state studied in the central field approximation where the 
electron motion is solution of the Dirac equation with an effective central potential 
created by a /c or a A; — x electron ion core; or any other atomic process involving 
highly excited electrons which needs to be studied using multichannel spectroscopy or 
quantum deffect theory (Owono Owono et al 2002, Mannervik et al 1998, Xiao-Min 
et al 1994, Aymar et al 1996). 

The Hamiltonians can be written as (Martinez-y-Romero 2000) 

=car[pr - i/3ei(ji + l/2)/r] + M^^c^ + Vi{r), 
H2 =car\pr - ipe2{j2 + l/2)/r] + M2pc^ + V2{r), 

in writing equation (1) we assume Mi ^ M2 since it is convenient for the calculations 
we perform for deriving the recurrence relations. The eigenstates of these radial 
Hamiltonians correspond to a definite value of the total angular momentum (orbital 
plus spin) J = L + S and of the quantum number e = (—1)-'^'^^/^, where, respectively, 
j is the total, I the orbital, and s the spin angular momentum quantum numbers, the 
operators 



are Dirac matrices in the two-dimensional representation appropriate for the radial 
problem we are discussing, and the subscript simply distinguishes between the systems 
pertaining to the different potentials. Each Dirac equation can be written as 



HkMr) = EkMr) (3) 

where the energy eigenvalues En^j^s^ = ^ind the corresponding eigenfunctions 
i^ukjkSki''^) = i^kii^) ^re assumed known. In this work we use atomic units (h = rUe = 
e = 1) so the masses in (1) are given in units of the electron mass. The radial 
eigenstates ipk{i^) are the purely radial spinors (Drake 1996, Greiner 1991) 
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= i(^-W); (4) 

F and G are, respectively, the big and the small components of the energy eigenstate. 
In the following we employ, as we did in equation(3), the shorthand k for describing the 
whole set of quantum numbers Ukjk^k- We also use the ket notation for the relativistic 
eigenstates Inkh^k) = 

2.1. The first recurrence relation for the case of two centres 

Taking the difference between the radial Hamiltonians Hi and in (1), we obtain 



Hi^H2 + icarP^ - ^(3M- - {V^ir) - Vi(r)) . (5) 
zr 

where = M2 ± Mi, and A=^ = €2(2^2 + 1) ± ei(2ji + 1). On employing (1) again, 
we can directly evaluate the commutator 



[HiJ{r)]^-icar^ (6) 

where /(r) is an arbitrary radial function and [iJ, /(r)] stands for the commutator 
between H and f{r). We can calculate this commutator again, but now using equation 
(5), to get the alternative form 

[i/i, /(r)] = H^f{r) - f{r)Hi + (^ica./3^ - c^/JM" - V'^ f{r). (7) 
It is now simple to obtain, from equations (6) and (7), the relation 



{E2 - Ei){2\f\l) = (2| {c'/3M- + V-) f\l) - ic{2\ar (^f + P^f^ (8) 

where we have additionally taken matrix elements between the eigenstates (4), we use 
the notation (1| = (nijiei| and |n2j2e2) = |2), and we have defined 



V^ = V2{r)±Vi{r). (9) 

Equation (8) may lead to recursions between relativistic matrix elements of radial 
functions between hydrogenic states (Martinez-y-Romero et al 2003, Ilarraza-Lomeli et 
al 2002, Martinez-y-Romero et al 2002, Martinez-y-Romero et al 2001) and generalizes 
a nonrelativistic one useful for similar purposes (Nunez- Yepez et al 1995). Equation 
(8) is an exact relation for the calculation of any /(r) two-centre matrix elements in 
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relativistic quantum mechanics. Taking the potentials as equal, i.e. Viir) — V2(r), we 
recover a relation which has proved useful for obtaining recurrence relations between 
atomic matrix elements in relativistic quantum mechanics (Ilarraza-Lomeli et al 2002, 
equation (6); Martinez-y-Romero et al 2001, equation (24)). 

Albeit exact, equation (8) is not entirely convenient due to the presence of the 
operator arP. To get rid of this factor, we found it convenient to deal directly with 
operator relations and not with the matrix elements themselves. The matrix elements 
will be evaluated at the end of the operator calculations. 

Let us first establish that 



H2f - fH, = {c'pM- + V-)f- icar (^f + pf^^ , 



(10) 



notice that equation (8) above can be obtained from (10) just by taking matrix 
elements. The following result is also easily established 

H2f + fH, = (c2/3M+ + y+) / - icar (2/-^ + f' + ^ + Pf^^ (H) 



Then, it can be seen that 



- ic {H^arf + ar fH,) = ica, (c^/JM" - V+) / - (^2/|- + /' + ^ - , (12) 



and that 



H2fV- - fV-H, = {c^(3M- + V-) V-f - icar ( V-f + ^/ + PfV-^ ) . (13) 



The next relation is also readily apparent 



ic 



/ / 
H2ar(3 — h ar(3—Hi 
r r 



= -icarpV+--c^ 



r r^ / 2r r 



■ (14) 



To further proceed, let us define 'ip{r) = H2f{r) + f{r)Hi, and evaluate 



iy^V - i^H^ = c'(3^f + c' (^^^ f + {c'(3M- + F") V 



2r \ ar r 



—icar 



f + {V- - c'pM^) + c'pM- (^2/|- + /' + ^) 



(15) 
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Given these last expressions [equations (10)-(15)], it is relatively simple to obtain, 
from (10), 



icar { f + Pf— ] = {H2f - fHi) - {c'pM- + V-) /, 



(16) 



and from (11), 



ICCKr 



rP (2f^ + f+K^= (H^Pf + PfH^) - {c'pM+ + V+) /3f + icar^f. (17) 
From equation (12) we obtain 



2r \ ar r 



A- A- 
- ic ( H2ar(3—f + ar(3—fHi 



-ica. 



2r 



(18) 



From equation (13) we obtain 



-icq;. 



{v-f + = {H2fV- - fV-H,) - {c'pM- + V-) V-f 



+icarV-—Pf. 
2r 

Substituting equations (16), (17), (18), and (19), into equation (15), we get 



-c'M+ {H^Pf - PfH,) + c^M+V-pf + c'M- {H^Pf + PfH,) 
-c'M-V+Pf + V-[2 {H2f - fHi) - V- ] 

-ic (H2arP^f + arP^fH^ - icar^ {c'pM- - V+) pf. 



(19) 



(20) 



The terms which include the operator icar can be also obtained from (14), thus we 
get the expression 



A A 

-ic ( H2arP—f + arP—fHi 



-icar— {c'PM- - V+) Pf 



—c 



2r 



P f- 



f 



2r 



■f 



(21) 
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Susbtituting now equation (21) into (20), it easily yields 



H,i; - ^H, = -c' (/" - Pf^^ + 2c'^^Pf + {M-fc'f 
-c^M^ {H2(3f - (3fH^) + c^M+V-(3f + c^M" {H^f^f + (3fHi) (22) 
-c'M-V+Pf + V-[2 {H,f - fH^) -V-]-^^ [(3f - . 

Evaluating the matrix elements between the Dirac eigenstates (2| and |1) and 
rearranging, we finally obtain the relation 

ao(2|/|l) + a2(2|^|l) - 2E-{2\V-f\l) + {2\{V-f f\l) + c^2\f"\l) 



= bo{2\pf\l) + b,{2\p^\l) - c'M-{2\V+Pf\l) + c'M+{2\V-pf\l) + b,{2\pL\i)^ 
where 



ao = {E')'-{c'M-Y 

a-2 = - -A-A+ 
4 

bo =S {M-E+ - M+E-) (24) 
bi =c'A- 

h =^ (A+ - A-) 

This is the first relation between matrix elements of an arbitrary radial function /(r) 
between eigenstates of two different potentials as a function of the eigenenergies in 
relativistic quantum mechanics. 

2.2. The second two-centre recurrence relation 

Given that the radial eigenstates have two components in relativistic quantum 
mechanics, it should be clear that we need more relations. To obtain such second 
equation, let us evaluate the following operator identity — i.e. again we calculate 
with the basic operators before any matrix element is taken. 

H2fV- + fV-H, = {c^pM+ + V+) V-f 

\ ar ar r 2r J 

Using again the definition 'ip{r) = H2f{r) — f{r)Hi, we obtain 
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+ {c'pM-^ + {c'pM- + V-)f- (^2/'^ + r + 2^ - (26) 

-icar i'^V-f^ + V-f + ^/ + 2F-^ + + c^PM-f + c^M-^/) . 

\ ar ar r Ir 2r J 

The calculations required for getting to the recurrence relation are similar to that 
used in the last section. So, from equation (10) we obtain 

- icarc'pM- (^f + (3f^^ = c^M- {H,f - fH,) - c^M- {c^pM' + F") Pf, (27) 
from (12) we get 



-icar {c'pM- - V+) /' - c'pf^, 
and from equation (25), we obtain 

-icar (^V-f^ + V-f + + 2V-^) = {H,fV- + fV-R,) 

\ ar ar r J 

+icarV-pf^ - {c^pM+ - V+) V-f. 
Using equations (16), (27), (28) and (29), in equation (26), it yields 

H.^ + ^H. = - c^p^ (/ - I) + c^^^f - V-V-f 

+M-M+c'f + V+ {H^f - fH^) + V- {H^f + fH^) 
-ic {H2arf' + arf'Hr) - icar {c^pM- - V+) /'. 

Again, the last two terms of (30) can be obtained from (14), and combining these with 
-ic {H2arf + arf'Hi) -icar {c^PM- - V+) f = 



(28) 



(29) 



(30) 
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we get 



+M-M+c''f + V+ {H2f - fH,) + V- (H^f + fH,) (32) 



r J 2r 

Taking matrix elements, we obtain the recurrence relation 

Co(2|/|l) + a2(2|^|l) - E+{2\V-f\l) - E-{2\V+f\l) + {2\V+V-f\l) 



r^ 

c'{2\^l) + c^2\ni) = |(2|/34|1) + h{2\P^\l), 



(33) 



where the only newly defined coefficient is 

Co = E^E- - c^M^M-. (34) 

2.3. The third and fourth two-centre recurrence relations 

To get the third recurrence relation, we have to substitute, in a similar way to 
what we have done for getting the previous two relations, equations (10), (12), (14) 
and (25) into (26) and, after some juggling with the resulting terms, and taking matrix 
elements, we finally obtain the third recurrence relation 

eo(2|/|l) = goWf\l) - (2| {V^ - V') /3/|l), (35) 

where 

Co = (M+ - M-) 

^ ' (36) 

g^ = E+ -E-. 

This is a very simple equation that, besides, allows writing the matrix elements of / 
in terms of those of [3f . To take advantage of this fact, substitute equation (35) into 
(23) to obtain a new relation 

ylo(2|/3/|l) + ^i(2|/3^|l) + A2{2\V-^f\l) + ^3(21 {V-f /3/|l) 
+ (2| {V-f^f\l) + A,{2\V^Pf\l) + 2E-{2\V-V^^f\l) 
- (2| {V-y V+Pf\l) + a2{2\ {V+ - V-) = 

A,{2\pi^\l) + c'go{2\pf"\l) - c'{2\ {V+ - V") pf"\l). 
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where the newly defined coefficients are 





{E-Y{E- 


Ai = 


-T - 


A2 = 


-2E- {E+ 


A,= 


E+ - 3E- 


A,^ 


c^M+M- - 


A,= 


-(M+-1 
2 ^ 



r-\2 



(38) 



Equation (37) is the fourth recurrence relation for the calculation of relativistic /(r) 
two centre matrix elements in terms of the energy eigenvalues of the intervening 
potentials. Notice that, at difference of the previous relations [equations (23), (33) 
and (35)], equation (37) just relates among themselves matrix elements of j3f and its 
derivatives times a certain function of r. 



3. Conclusions 



We have obtained recurrence relations for the calculation of two-centre matrix elements 
of a radial function between states of two different radial potentials. The obtained 
recursions are given in the most general case of an arbitrary function taken between any 
non necessarily diagonal radial eigenstates of two radial potentials. These recursion 
relations have, as particular cases, recursions between one-centre integrals or, in other 
particular cases, between overlap and one centre integrals in Dirac relativistic quantum 
mechanics. We expect the obtained recursions, together with the previous one-centre 
relations we have obtained (Ilarraza-Lomeli et al 2002, Martinez- y- Romero et al 2001, 
2002), to be useful in atomic or molecular physics calculations as they may simplify 
calculation in the range of applicability of Dirac's relativistic quantum mechanics 
(Bang and Hansteen 2002, Moss 1972). For most uses of the relations we first have to 
set Ml = M2, i.e. M~ = and M+ = 2 — if the particles arc electrons; since the use 
of unequal masses is just a recourse of our calculational method. 

From a practical angle, there is little that can be done for the analytical evaluation 
of two centre integrals of atomic physics interest beyond the Coulomb and the few 
similarly exactly solvable potentials. However, there are numerical methods that, 
after being adapted to relativistic conditions, can provide the crucial "seed" results 
needed for the systematic use of the recurrence relations obtained here (Chen et 
al 1993). Our results can be also useful in the so-called perturbation theory 
of relativistic corrections, in relativistic quantum deffect calculations, and for the 
relativistic extension of the calculations of exchange integrals using Slater orbitals or 
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Coulomb-Dirac wave functions (Owono Owono et al 2002, Bang and Hasteen 2002, 
Charro et al 2001, Kuang and Lin 1996, Rutkowski 1996). It is also possible that our 
relations could be applied to a generalization of an approximate two-centre technique 
used for studying electron-impact ionization of simple molecules (Week et al 2001). 

On the other hand, our results can be also of interest in nuclear studies. Since the 
3D Woods-Saxon potential, used in the Dirac equation for describing the interaction of 
a nucleon with a heavy nucleus, has been explicitly solved recently and its eigenfuctions 
expressed in terms of hypergeometric functions (Jian-You et al 2002), so it has the 
features needed for the direct use of our recurrence relations. 
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